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Abstract 

An approach to deep inelastic scattering is described in which the matrix elements arising 
from the operator product expansion are factorised into composite operator propagators 
and proper vertex functions. In the case of polarised \ip scattering, the composite operator 
propagator is identified with the square root of the QCD topological susceptibility y/xUfi), 
while the corresponding proper vertex is a renormalisation group invariant. We estimate 
x'(0) using QCD spectral sum rules and find that it is significantly suppressed relative to 
the OZI expectation. Assuming OZI is a good approximation for the proper vertex, our 
predictions, ft dxg{{x;Q 2 = lOGeV 2 ) = 0.143 ±0.005 and G^ 0) = AE = 0.353 ±0.052, are 
in excellent agreement with the new SMC data. This result, together with one confirming 
the validity of the OZI rule in the rj' radiative decay, supports our earlier conjecture that 
the suppression in the flavour singlet component of the first moment of g\ observed by 
the EMC-SMC collaboration is a target-independent feature of QCD related to the U(l) 
anomaly and is not a property of the proton structure. As a corollary, we extract the 
magnitude of higher twist effects from the neutron and Bjorken sum rules. 
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1. Introduction 



The discovery by the EMC collaboration [1] (see also [2]) of an unexpected suppression 
of the first moment of the polarised proton structure function g\ has provoked an extensive 
discussion of the parton model interpretation of QCD in deep inelastic scattering processes 
involving the axial U(l) anomaly. (For reviews, see refs.[3,4].) While it has so far proved 
possible with careful redefinitions and interpretations [5] to preserve the essence of the 
parton model description, it is becoming clear that these processes involve subtle field 
theoretic properties of QCD which lead beyond both the original and QCD-improved 
parton approximation. In this paper, we develop an alternative approach to deep inelastic 
scattering emphasising field theoretic concepts such as the operator product expansion 
(OPE), composite operator Green functions and proper vertices. This clarifies some of the 
difficulties encountered in the parton description and gives a new insight into the underlying 
reason for the EMC result. In particular, our analysis strongly suggests that the observed 
suppression of the first moment of g\ is a generic QCD effect related to the anomaly and is 
actually independent of the target. Rather than revealing a special property of the proton 
structure, the EMC result reflects an anomalously small value of the first moment of the 
QCD topological susceptibility [6, 7]. 

The essential features of this method are easily described for a general deep inelastic 
scattering process. The hadronic part of the scattering amplitude is given by the imaginary 
part of the two-current matrix element {N\ J M (<z) J v {—q)\N) illustrated in Fig. 1, where 
is the current coupling to the exchanged hard photon (or electroweak vector boson) and 
\N) denotes the target. The OPE is used to expand the large Q 2 limit of the product of 
currents as a sum of Wilson coefficients Ci(Q 2 ) times renormalised composite operators 
Oi as follows (suppressing Lorentz indices), 



The dominant contributions to the amplitude arise from the operators Oi of lowest twist. 
Within this set of lowest twist operators, those of spin n contribute to the nth moment of 
the structure functions, i.e. 



The Wilson coefficients are calculable in QCD perturbation theory, so the problem 
reduces to evaluating the target matrix elements of the corresponding operators. We now 
introduce appropriately defined proper vertices Fqnn, wmcri are chosen to be 1PI with 
respect to a physically motivated basis set Ok of renormalised composite operators. The 
matrix elements are then decomposed into products of these vertices with zero-momentum 
composite operator propagators as follows, 
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(1.2) 




(1.3) 
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This is illustrated in Fig. 2. In essence, what we have done is to split the whole ampli- 
tude into the product of a "hot" (high momentum) part described by QCD perturbation 
theory, a "cold" part described by a (non-perturbative) composite operator propagator 
and finally a target-dependent proper vertex. 

All the target dependence is contained in the vertex function T^^^. However, these 
are not unique - they depend on the choice of the basis Ok of composite operators. This 
choice is made on physical grounds based on the relevant degrees of freedom, the aim 
being to parametrise the amplitude in terms of a minimal, but sufficient, set of vertex 
functions. A good choice can often lead to an almost direct correspondence between the 
proper vertices and physical couplings such as, e.g., the pion-nucleon coupling g^NN- In 
particular, it will be wise to use, whenever possible, RG-invariant proper vertices. 

Despite being non-perturbative, we can frequently evaluate the composite operator 
Green functions using a combination of exact Ward identities and dynamical approxima- 
tions (see sects. 2 and 3). On the other hand, because of the target dependence, the proper 
vertices are not readily calculable from first principles in QCD, so we are in general left 
with a parametrisation of the amplitude in terms of a (hopefully small) set of unknown 
vertices. These play the role of the non-perturbative (i.e. primordial or not-yet-evolved) 
parton distributions in the usual treatment. Just as for parton distributions, many dif- 
ferent QCD processes can be related through parametrisation with the same set of vertex 
functions. 

Now compare this approach with the parton model. In the original parton model, 
the amplitude is approximated by Fig. 3, describing the scattering of a large Q 2 photon 
with a parton in the target nucleon. This picture is already sufficient to reveal Bjorken 
scaling. It may be improved in the context of QCD by including gluonic corrections, 
exactly as in the OPE, as shown in Fig. 4. These give the logarithmic scaling violations 
characteristic of perturbative QCD. The total amplitude is therefore factorised into a 
perturbative scattering amplitude for the hard photon with a parton (quark or gluon) and 
a parton distribution function giving the probability of finding a particular parton with 
given fraction x of the target momentum. 

The question of whether the full QCD amplitude can be given a natural parton inter- 
pretation depends on the composite operators Oi in the Wilson expansion. For example, 
if the lowest twist operator for a given process is multilinear in the elementary quark and 
gluon fields rather than simply quadratic then the diagram of Fig. 4 is not appropriate and 
the process can only be described in terms of multi-parton distributions [8]. A more subtle 
problem arises when the operators Oi are non-trivially renormalised and mix with other 
composite operators under renormalisation. In this case, the parton interpretation is pre- 
served by defining parton distributions directly in terms of the operator matrix elements 
(see, e.g., ref.[8]). This procedure becomes especially delicate [5] in the case of polarised 
deep inelastic scattering because of the special renormalisation properties of the relevant 
Wilson operator J® 5R due to the axial U(l) anomaly. 

In this paper, rather than attempt to interpret the amplitude for polarised deep in- 
elastic scattering in terms of specially defined polarised quark and gluon distributions, we 
instead focus the analysis on the composite operator level. By splitting the matrix elements 
in the form of eq.(1.3), we can exploit chiral Ward identities and the renormalisation group 
to separate out generic features of QCD manifested in the composite operator propagator 
from specific properties of the target. In the next section, we see how this clarifies the 
origin of the suppression of the first moment of g\ observed in polarised \xp scattering. 
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2. The First Moment Sum Rule for g\ 

Our starting point is the familiar Ellis-Jaffe[9] sum rule for the first moment of the 
polarised proton structure function g\. For Np = 3 and in the MS scheme [10], this reads*: 

T\(Q 2 ) = f dxg{{x;Q 2 ) 
Jo 

(0) + ^ 8) (0)) (l - f - 3-583(^) 2 - 20.215(^) ! 



+ ^(o ; g 2) (i-^-o,5o( ? ) 



(2.1) 



where the are form factors in the proton matrix elements of the axial current 



(P\j; 5R (k)\P) = G { 2 } (k 2 )u W u + G^\k 2 )k^u, (2.2) 
and a is an SU(3) flavour index. In our normalisations (see ref.[7]): 

Gf = \{Au - Ad) 

G { * } = -^=(Au + Ad- 2 As) 
2y 3 

G ( ° ] = Au + Ad + As = AS. (2.3) 

We ignore heavy quarks and, for simplicity, set the light quark masses to zero in the 
formulae below. 

The axial current occurs here since it is the lowest-twist, lowest spin, odd-parity 
operator in the OPE of two electromagnetic currents, i.e. 

JMM-q) ~ 2 J2 ^^C a (Q 2 ) J a p5R + .... (2.4) 

Q a =0,3,8 V 

The suffix R emphasises that the current is the renormalised composite operator. Under 
renormalisation, the gluon topological density Qr and the divergence of the flavour singlet 
axial current J® 5R mix as follows[12], 

Qr = Qb - ^ 1 - Z ) d " J %B, (2-5) 



* We use the NLO and NNLO coefficients given in Ref. [11]. However, due to our definition (2.5) 
of the renormalised composite operators, the radiative corrections of the singlet are different from the 
corresponding terms in ref. [11], which uses a different renormalisation of the singlet operators. 
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where J® 5B = YIqIhIsI an d Qb = f^trG^G^ and we have quoted the formulae for N F 
flavours. The mixing is such that the combination occurring in the axial anomaly Ward 
identities, e.g. 

(0\(d»J% R -2N F Q R ) 6 k \0) + (0\5 5 6 k \0) = 0, (2.6) 
is not renormalised. 

Since J® 5R is renormalised, its matrix elements satisfy renormalisation group equations 

with an anomalous dimension 7, so that in particular G^ (0; Q 2 ) depends on the RG scale 
(which is set to Q 2 in eq.(2.1)). 

As we have emphasised elsewhere, G^ does not, as was initially supposed, measure the 
spin of the quark constituents of the proton. The RG non-invariance of J® 5R (a consequence 
of the anomaly) is itself sufficient to prevent this identification. The interest in the first 

EMC data [1,2] on polarised fxp scattering*, which allows the following result for G^ to 
be deduced, 

G { ^(0;Q 2 = HGeV 2 ) = AE = 0.19 ± 0.17, (2.7) 

is rather that this value for G^ represents a substantial violation of the OZI rule[13,14], 
according to which we would expect 

G { a\0)ozi = 3F-D~ 0.579 ± 0.021. (2.8) 
Here, we have used [15,16]: 

F + D ~ 1.257 ±0.008 F/D ~ 0.575 ± 0.016 (2.9) 



as fitted from hyperon and (3 decays. The assumption that the OZI rule is satisfied for 
Ga(Q) is equivalent to the Ellis-Jaffe sum rule prediction for the first moment of g\. 

It follows immediately from eq.(2.2) (assuming the absence of a massless pseudoscalar 
boson in the U(l) channel) that 

4 0) (0;Q 2 ) u ToU = JL (p\&>j% r \p), (2.10) 

where M is the proton mass. The anomalous chiral Ward identity then allows G^ to be 
re-expressed as the forward matrix element of the renormalised gluon topological density 
Qr, i-e. 

Gf{^Q 2 ) u l5 u = ^-2N F (P\Q R (0)\P). (2.11) 



* The combined SLAC/EMC data quoted in ref.[l] gives 

rf(Q 2 = HGeV 2 ) = 0.126 ±0.010 ±0.015 

The result for in eq.(2.7) is extracted from the sum rule using the values for F and D given above 

and the running coupling from tau-decay data[30] (see the remarks after eq.(3.32)). 



4 



Notice that in terms of bare fields, Qr contains both gluon and quark bilinears. This, 
together with the explicit factor of a s in the definition of the topological density, is the 
source of the difficulty in giving a natural and unambiguous parton interpretation^, 3, 4]. 

At this point, we apply the method described in the introduction. We choose as the 
composite operator basis Ok the set of renormalised flavour singlet pseudoscalar operators, 
viz. Qr and and Q5R, where, up to a crucial normalisation factor discussed below, the 
corresponding bare operator is simply the singlet i (Z75(Z- We then define T[Qr, $5i?; P, P] 
to be the generating functional of proper vertices which are 1PI with respect to these 
composite fields only. (Here, P and P denote interpolating fields for the proton - they 
play a purely passive role in the construction.) V is obtained from the QCD generating 
functional by a Legendre transform with respect to the sources for the composite operators 
Qr and $5i? only. We may then write (c.f. eq.(1.3)) 

(p\q r (o)\p) = (o\Qr(o) q r (o)\o) r QRP p + (o\Q R (o) *5*(o)|o> r„ 5KP p, (2.12) 

where the propagators are at zero momentum. 

The composite operator propagator in the first term in eq.(2.12) is the zero-momentum 
limit of an important quantity in QCD known as the topological susceptibility x(/c 2 ), viz. 

X(k 2 ) = Jdxe ik - x i(0\T*Q R (x)Q R (0)\0). (2.13) 

The second term is clearly independent of the normalisation of the renormalised quark 
bilinear operator $5i?. We choose to normalise this operator in such a way that the inverse 
two-point function r$ 5W $ 5H , which has to vanish at k 2 = 0, is equal to k 2 , the correct 
normalisation for a free, massless particle. With this normalisation, a straightforward 
but intricate argument [7] using chiral Ward identities (see Appendix A) shows that the 
propagator (0\Qr &sr\0) at zero momentum is simply the square root of the first moment 
of the topological susceptibility x(/c 2 ). We therefore find, 

(P\ Q R (0) |P> = X (0) T QrPP + v^(0) (2.14) 

The chiral Ward identities further show that for QCD with massless quarks, %(0) ac- 
tually vanishes. (This is in contrast to pure Yang-Mills theory, where x(0) is non-zero and 
is related to the rj' mass in the large Nc resolution of the U(l) problem[17,18].) Only the 
second term in eq.(2.14) remains. Remarkably, this means that the matrix element of the 
renormalised gluon density Qr measures the coupling of the proton to the renormalised 
pseudoscalar quark operator $5^. This happens because the composite operator propaga- 
tor matrix in the pseudoscalar (Qr, $>5r) sector is off-diagonal. We therefore arrive at our 
basic result [7] , 

G^(0;Q 2 )u l5 u = ^2N F ^W)Y^ rPP . (2.15) 

The renormalisation group properties of eq.(2.15) are central to our argument. With 
the normalisation of chosen above, it can be shown[7] that the proper vertex T§ P p 
is RG invariant and so has no scale dependence. The scale dependence needed to match 
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Gjp is provided entirely by the topological susceptibility which, as shown in Appendix A, 
satisfies the RGE 

("! + « a >«^-K)*'( o > = °- < 2i6 » 

The challenge posed by the EMC data is to understand the origin of the OZI violation 
in Gjl - The OZI approximation applied to the r.h.s. of eq.(2.15) would require* (neglecting 
flavour 577(3) breaking) T^ BRP p ~ V^g m NN u^ 5 u while a/x'(0) ~ ( 1/^/6)/^ • 

Our proposal is that we should expect the source of the OZI violation to lie in RG 
non-invariant terms, i.e. in %'(()). The reasoning is straightforward. In the absence of the 
U(l) anomaly, the OZI rule would be an exact property of QCD. So the OZI violation is 
a consequence of the anomaly. But it is the existence of the anomaly that is responsible 
for the non-conservation and hence non-trivial renormalisation of the axial current </° 5Ji - 
We therefore expect to find OZI violations in quantities sensitive to the anomaly, which 
we identify through their RG dependence on the anomalous dimension 7. This seems 
reasonable since, if the OZI rule were to be good for such quantities, it would mean 
approximating a RG non-invariant, scale-dependent quantity by a scale-independent one. 



If this proposal is correct, we expect a/x'(0) to be significantly suppressed relative to 
its OZI approximation of (l/v^})/^. The proper vertex T§ PP would behave exactly as 
expected according to the OZI rule. That is, the Ellis-Jaffe violating suppression of the first 
moment of g\ observed by EMC would not be a property of the proton at all, but would 
simply be due to an anomalously small value of the first moment of the QCD topological 
susceptibility %'(0). 

In the next section, we attempt to verify this hypothesis by evaluating x'(0) using 
QCD spectral sum rules. 



3. QCD Spectral Sum Rule Estimate of x'(0) 

We now present an estimate of x'(0) in QCD with massless quarks using the method 
of QCD spectral sum rules (QSSR) pioneered by Shifman, Vainshtein and Zakharov [19] 
and reviewed recently in ref.[20]. 

The correlation function x(k 2 ) is defined in eq.(2.13) and its renormalisation group 

* To understand this, we note from ref.[7] that eq.(2.15) is equivalently written as one form of the U{1) 
Goldberger-Treiman relation, viz. 

G^^Q 2 ) = F VOZI g VOZI NN 

where F rjQZI and g VOZI NN are respectively the decay constant and nucleon coupling of a state \rjozi}- 
\VOZi) is an unphysical state in QCD (i.e. not a mass eigenstate) which in the OZI or large Nc limit, 
in which the anomaly is absent, can be identified as the massless U(l) Goldstone boson. Simple quark 
counting rules then relate g-q OZI NN to the 778-nucleon coupling g Vs NN- 

This identification is the origin of our choice of normalisation of &sr. In the OZI limit, r$ 5J? jvJv 

becomes the Goldstone boson - nucleon coupling. 
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equation is given in Appendix A. Including the inhomogeneous contact term[21], we have 



with the beta function 



a. q djJL V 7T / V 7T / 



(3.2) 



where, for QCD with N F flavours, (3 X = -5 (11 - and /3 2 = - J (51 

the anomalous dimension[12] 



f iV F ), and 



7 = fi-^-logZ = -(—) 

a/j, V 7T / 



(3.3) 



The extra RG function /3^ L ) (so called because it appears in the longitudinal part of the 
Green function of two axial currents) is given by 



1 /?W 
(2iV F )2 P 



1 /« s \ 2 



327T 2 V 7T 

The RGE is solved in the standard way, giving 

-2 f*dtr>r(a a (fj) 



1 + 



29 /a 



4 V 7T 



(3.4) 



x(/c ,a a ;/i) 



x(/c 2 ,a s (t);^e*) 
- 2 f dt" ^ L \a s {t")) e 2 /o'W«-<*'>) 



, (3.5) 



where ct s (t) is the running coupling. 



The perturbative expression for the two-point correlation function in the MS scheme 



is[22] 



X(k 2 ) 



P.T. 



/a s \ 2 2 4 -k 2 



1 + 



log 



-A; 2 29 



+ — +... 



(3.6) 



The non-perturbative contribution from the gluon condensates (coming from the next 
lowest dimension operators in the OPE) is [23] 



a. 



16tt 2 



1 _ / CK 



Z V 7T 



log 



/// 



O:, 



(3.7) 



The RGE has been used to check the consistency of the leading log approximation in 
the perturbative expression and to fix the radiative correction in the gluon condensate 
contribution. 
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For the QSSR analysis of y/(0), we use the subtracted dispersion relations 

dt 1 1 



f (*(*'> -*(»>) = I 



t t — k 2 — ie it 



Im X (t) 



and 



¥ (x(k 2 )-x(o)-k\'(o)) = j o 



dt 



1 1 



t 2 t — k 2 — ieir 



Imx(t). 



(3.8) 



(3.9) 



Then, taking the inverse Laplace transform[20] of both sides of the dispersion relations 
and using the fact that x(0) = in massless QCD, we find* 



I 



tc dt 1 

^e-^-lm X (t) 

t 7T 

-«,\2 2 



\8nJ 



,P2 



* ^r- 2 (l-e-*-(l+t c r)) 1 + (g) (^29 + 4/3^1 - 7s ) - 8^ log(- log rA 2 ) 



(3.10) 



and 



X'(0) * / 
Jo 



tc dt 



t 2 

- (s) 2 ^- 1 f 1 - I 1 + (s) ( 29 - ^ - 8 i log (- 1o ^ a2 » 



e tT — Imx(t) 



7T 



.ft 



+ 



(3.11) 



where ct s is the running coupling expressed in terms of the QCD scale A from the two-loop 
relation: 



TV 



7T 



i-(f)|M-iog^) 



with 



a. 



7T /?i log TA 2 ' 



(3.12) 



(3.13) 



In these expressions, we have cut off the t integration at some scale t c and used the 
perturbation theory approximation to Imx(t) for t > t c . 

In order to extract a value for x'(0) from these sum rules, we keep only the lowest 
resonance (the rj') contribution to the spectral function, i.e. we assume 



1 



n 



Im X (t) = 2mJ, mj,) + "QCD continuum" 9(t-t c ), (3.14) 



For the corresponding results in pure Yang-Mills theory, see refs. [24,25]. 
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where m v > is the mass of the rj' extrapolated for massless QCD, viz. 

m\, ~ m%, - \m K 2 ~ (0.87GeV) 2 . (3.15) 

To evaluate eqs.(3.10) and (3.11), we use 

A ~ (350±100)MeV (3.16) 

for the QCD scale parameter [26], 

(a s G 2 ) ~ (0.06±0.02)GeV 4 (3.17) 

from a global fit of the light mesons and charmonium data[20], and parametrise the triple 
gluon condensate as 

(g 3 G 3 ) ~ (1.5±0.5)GeV 2 (a s G 2 ) (3.18) 

using the dilute gas instanton model[19]. We show the result in Fig. 5a for x'(0) plotted 
versus r for different values of t c . In Fig. 5b we show the behaviour of the r minima for 
different t c . Our optimal result corresponds to the range of values of t c corresponding to 
the first appearance of the r minimum until the beginning of the t c stability region. The 
value of r at which the stability occurs is around 0.4 to 0.6GeV , which is quite small 
compared with the light meson systems and is consistent with qualitative expectations [23] 
of a scale hierarchy in the QSSR analysis of gluonium systems. This small value of r also 
ensures that higher-dimension operators such as those arising from instanton-like effects 
will not contribute in the OPE. We deduce: 

vV(0) ^ (22.3 ± 3.2 ± 2.8 ±1.3)MeV, (3.19) 

where the first error comes from (a s G 2 ), the second one from A and the third from the range 
of t c values from 4.5GeV 2 to 7.5GeV 2 . The effects of the triple gluon condensate and the 
radiative corrections are relatively unimportant, contributing about 3 — 10% to x'(0). We 
add a guessed error of 5% each from the unknown non-perturbative and radiative correction 
terms. Finally, adding all these errors quadratically, we find the following Laplace sum rule 
estimate of the first moment of the topological susceptibility evaluated at r ~ 0.5GeV -2 : 

vV(0) ^ (22.3±4.8)MeV. (3.20) 

As a check on the validity of this result, we now repeat the analysis using the finite 
enrgy sum rule (FESR) local duality version of the spectral sum rules discussed in ref.[25]. 
The advantage of the FESR method is that it projects out the effects of the operators of a 
given dimension[27] (in this case, dimension 4) in such a way that, at the order to which 
we are working, the FESR analogues of the sum rules (3.10) and (3.11) are not affected 
by higher dimension operators such as those induced by instanton-like effects. 
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The FESR sum rules are 

rtc dt 1 



J 

Jo 



2 
_c 

2 2 



,Pi 



+ 



V8tt 



) 



2tt 



(3.21) 



and 



X'(0) 



tc dt 1 T , , 

87T 7 7T 2 C 



1+ (£)( 29 - 4 ^ 



8^1og(-logrA 2 ) 
Pi 



(3.22) 



Analysing eqs.(3.21) and (3.22), we realise that the solution increases monotonically with 
t c so that no firm prediction can be made, although the result gives a rough indication 
of consistency with the previous Laplace sum rule. To overcome this problem, we repeat 
the analysis using only the FESR (3.22) and using as an extra input the value of the 
parameter f v / extracted from the first Laplace sum rule (3.10). The value of f v > is given 
in Appendix B. This weakens the t c dependence of the result and t c stability now appears 
as an inflection point. We obtain the result shown in Fig. 6 for different values of /„/ and 
A, from which we deduce that with t c ~ 6.5 — 9.5GeV , 



vV(0) ^ (25.5 ± 1.5 ± 2.0 ±1.0)MeV, 



(3.23) 



where the errors come from f v > , A and t c respectively. Adding the errors quadratically and 
including a further 5% error from the unknown higher order terms, we obtain at the scale 
t ~ 0.5GeV" 2 : 

vV(0) - (26.5±3.1)MeV, 

where, we have run the result from t c 
solution expressed in terms of A, viz. 



(3.24) 

8GeV 2 to the scale r = 0.5GeV" 2 using the RGE 



X'(0;^) ~ x'(0)e^, 



(3.25) 



where x'(0) is RG invariant. (Notice that the inhomogeneous term proportional to (3^ 
does not contribute to the first moment at k 2 = 0.) We see that the FESR result is 
consistent with the Laplace one. 

Taking the average of the Laplace and FESR results, we obtain our final estimate of 
the first moment of the topological susceptibility at the scale r = 0.5GeV -2 : 



vV(0) ^ (25.3±2.6)MeV. 



(3.26) 



This result should be compared with that obtained[24,25] in pure Nc = 3 Yang-Mills 
theory using a similar QSSR approach: 



^T7W)\ YM * (7±3)MeV. 



(3.27) 
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It is important to notice that this pure Yang-Mills result has been confirmed by lattice 
calculations[28,29], which is a strong indication of the validity of the methods used in 
deriving both (3.27) and (3.26). The introduction of massless quarks has changed the sign 
of x'(0) an d increased its absolute value by a factor of around 4. From the QSSR analysis, 
this effect is due mainly to the low value of the r/ mass of 0.87GeV (massless QCD) which 
enters into the spectral function, compared with the pseudoscalar gluonium mass of about 
1.36 - 1.66GeV [24,20] in pure Yang- Mills theory. 

To compare with the experimental result on the polarised proton structure function, 
we use the RGE to run the result for x'(0) to the EMC scale of lOGeV 2 . We find 

V7W)\ EMC * (23.2±2.4)MeV. (3.28) 

This is smaller by a factor of 1.64 ± 0.17 than the OZI value of (l/y/6)f n . We therefore 
do indeed find a significant suppression of %'(0) relative to its OZI value. 

To convert this result into a prediction for the singlet form factor, we take our fun- 
damental expression (2.15) for G^ and equate the proper vertex P p with its OZI 
expression given by the Goldstone boson-nucleon coupling. In this way, we obtain: 

Gf(0) = Gf(0)oz, ^g- (3.29) 

Using the value of G^ (0)ozi m eq.(2.8) and including an additional error of approx. 
10% for the use of the OZI approximation for the proper vertex, we arrive at our final 
prediction: 

G^ 0) (0;Q 2 = lOGeV 2 ) ~ 0.353 ± 0.052. (3.30) 
Substituting this result* together with 

G f = \(F + D), (3.31) 

into the first moment sum rule (2.1), using the values of F and D from eq.(2.9), and 
neglecting the higher twist terms (which are certainly negligible at Q 2 = 10 GeV 2 ) we 
deduce: 

T^lOGeV 2 ) ~ 0.143 ± 0.005. (3.32) 

Here, we have used the coupling a s (m T ) = 0.347±0.030 extracted from tau-decay data[30]. 
One should notice that the radiative corrections decrease the leading order result by about 
12%. 



In terms of the quantities A«, Ad and As defined in eq.(2.3), we have at Q 2 = lOGeV 2 : 
Am = 0.84 ±0.01 Ad = -0.41 ± 0.01 As = -0.08 ± 0.02 
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Our result, eqs.(3.30) and (3.32), certainly goes in the right direction, i.e. that of 
reducing the prediction from the OZI (Ellis-Jaffe) value. At the time we obtained it, 
however, eq.(3.30) still appeared too high compared to the experimental result (2.7), which 
would have implied further OZI violations in the proper vertex. Amusingly enough, while 
this paper was being completed we learned of the new results from the SMC collaboration 
which, combined with the earlier proton data, gives the new world average [31]: 

rf (10 GeV 2 ) = 0.145 ± 0.008 ± 0.011, (3.33) 

from which we deduce 

G^ 0) (0; Q 2 = 10GeV 2 ) = AE = 0.37 ± 0.07 ± 0.10. (3.34) 
These results are now in excellent agreement with our predictions. 



4. Tests of the Bjorken sum rule and estimate of higher twist effects 

Recently, the SMC collaboration at CERN[31],[32] and the E142 collaboration at 
SLAC[33] have produced data on the polarised neutron structure function g™. Since our 
proposal requires that the flavour singlet suppression is identical for the proton and neu- 
tron, we see no reason why the Bjorken sum rule [34]: 

6T ^ = r p - r? = J dx q 2 ) - (* ; q 2 )) 

- \* - a i - ^ivT - ^ <«) 

should not hold, at least up to flavour SU(2) breaking. Provided the measurements are 
at sufficiently high Q 2 , the higher twist corrections related to the coefficients a p — a n can 
be neglected. Analysis of the combined proton and deuteron data as performed in ref.[35] 
gives at Q 2 = 5GeV 2 [31]: 

5Y p ~ n ~ 0.203 ±0.029, (4.2) 
to be compared with the QCD prediction, with a s (5GeV 2 ) = 0.32 ± 0.02, of 

5Y p ~ n ~ (0.176 ± 0.003) + ap ~" n , (4.3) 
From this, one can deduce the difference of the higher twist coefficients (in units of GeV 2 ): 



a p -a n ~ 0.135 ±0.145. 
12 



(4.4) 



We can pursue an analogous analysis for the first moment of the neutron structure 
function, which satisfies the sum rule (c.f. eq.(2.1)): 



7(Q 2 ) = f dxg?(x;Q 2 ) 
Jo 

Y-G<?>(0) + ^Gf(0)\ (l - ^ - 3-583(^) 2 - 20.21 8 (^) : 



+ (4-5) 

where we have included the higher twist contribution. Evaluating this quantity at Q 2 = 
2GeV 2 , where the SLAC data are available, we find 

I7(2GeV 2 ) ~ -(0.031 ± 0.006) + (4.6) 

Comparing this with the SLAC data[33], 

I7(2GeV 2 ) ~ -(0.022 ± 0.011), (4.7) 

and using eq.(4.4), we can extract the coefficients of the higher twist terms. In units of 
GeV 2 , we find: 

a p ~ -0.117 ±0.145 

a n ~ 0.018 ±0.025. (4.8) 

These values of the higher twist terms are consistent with the previous determina- 
tions[36,37] from QCD spectral sum rules. However, these sum rules would be affected 
by a more general choice of the nucleon interpolating field[20] (the one used in refs. [36,37] 
is not the optimal one) and by the well-known[20] large violation by a factor 2-3 of the 
vacuum saturation of the four-quark condensate, which is assumed in refs. [36,37] to be sat- 
isfied to within 10-20%. In addition, radiative corrections, which are known to be large in 
the baryon sum rules[20], can also be important here. More accurate data on the Bjorken 
and neutron sum rules, and/or a measurement of the proton sum rule at lower Q 2 , are 
needed to improve the results in eq.(4.8), which are necessary to test the validity of the 
QCD spectral sum rule predictions in ref.[37]. 
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5. Further Discussion 

In this paper, we have presented evidence that the experimentally observed suppres- 
sion of the first moment of the polarised proton structure function g\ (the so-called EMC 
"proton spin" crisis) is a target-independent effect reflecting a suppression of the first mo- 
ment of the QCD topological susceptibility x'(0) relative to the OZI expectation. Not only 
does (0) not measure the quark spin - its suppression is not even a property of the 
proton structure. 

It would be interesting to test this hypothesis directly by polarised deep inelastic scat- 
tering experiments on other targets not simply related to the proton by flavour symmetry. 
We have already studied the case of a photon target and have presented elsewhere [38] a 
new sum rule for the first moment of the polarised photon structure function gj measur- 
able in polarised e + e~ colliders. However, this turns out to be a special case because the 
electromagnetic U(l) anomaly contributes at leading order and so the gj sum rule does 
not display the suppression mechanism described here. Another possibility is to consider 
semi-inclusive processes in which a particular hadron with a fraction z of the incoming mo- 
mentum is observed in the target fragmentation region. It was recently suggested [39] that 
such cross-sections should be described in terms of new, non-perturbative hybrid functions 
M{z, x, Q), called "fracture functions". To the extent that an OPE can be used, it would 
be possible to represent M in terms of the forward matrix element of a composite operator 
between a suitable proton-plus-hadron state. In this case, one would again factorise M 
into a composite propagator of the usual type and a proper vertex involving four external 
hadron legs. If the suppression of the polarised structure function indeed originates from 
the propagator, as we suggest, such a suppression should also be found at the level of the 
(less inclusive) fracture functions. 

So far, we have only considered the first moment of g\. Of course, we would like to 
extend our approach to higher moments and discuss the full x-dependence of the structure 
function. This would require knowledge of the renormalisation properties and composite 
operator Green functions of the higher-spin axial currents and gluon densities [40], together 
with the associated proper vertices. 

Another possible line of development would be to try to develop techniques to estimate 
the proper vertices themselves, rather than just the composite operator Green functions. 
To the extent that the quenched approximation may be trusted for the proper vertices, 
lattice calculations could already be suitable for the task, and QCD spectral sum rule 
techniques could be used in conjunction to check the validity of that approximation. We 
recall that, in contrast, use of the quenched approximation directly for the matrix elements 
of the operator Q can be shown to be completely unreliable since these are affected by low- 
lying poles that should disappear after dynamical quark loops are added. This is another 
example of how the apparent complication introduced by our splitting of matrix elements 
into propagators and proper vertices can ultimately pay dividends. 

Finally, it would be interesting to attempt to apply this analysis of deep inelastic scat- 
tering using proper vertices to other QCD processes normally described in the language of 
the parton model rather than in terms of the OPE. Semi-inclusive deep inelastic scattering 
is one such example, but many other interesting possibilities can be considered, especially 
in the context of hadron-hadron collisions. 
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Appendix A Chiral Ward Identities and the Renormalisation Group 



The anomalous chiral Ward identities for Green functions of the pseudoscalar opera- 
tors Q R and are (for zero quark masses) 

ik„ (0\J° 5R (k) Q R (-k)\0) - 2N F (Q\Q R (k) Q R (-k)\0) = (Al) 

ikp (0\J° 5R (k) $ 5R (-k)\0) - 2N F (0\Q R (k) * 5R (-k)\0) + (0\6 5 $ 5R (-k)\0) = 0. 

(A.2) 

So, at zero momentum, assuming there is no physical massless U(l) boson, 

(0\Q R (0) Q R (0)\0) = 0, (A3) 
showing that the topological susceptibility %(0) vanishes for massless QCD, and 

<0|Qh(0) $ 5fl (0)|0> = -^L-2(* R ), (A4) 

where (& R ) is the VEV of the scalar partner of and is non-vanishing because of the 
quark condensate. 

The field &s R is normalised such that the two-point proper vertex T^ 5R ^ 5R = k 2 . 
This means that r$ 5H $ 5H is (minus) a component of the inverse propagator matrix in the 
pseudoscalar sector, i.e. 

r* B ** 8 * = (0\Qr Qr\0)({0\Qr$ 5 r\0) 2 - (0\QrQr\0){0\®5R®5r\0)) ■ (A5) 



Expanding to lowest order in k 2 gives 

r* 5R * 5 « = X(0) (0\Q R (0) $5i?(0)|0)- 2 k 2 + 0(fc 4 ), (A6) 
where we have written (0\Q R (k) Q R (-k)\0) = x'(0)/c 2 + 0(k 4 ). We therefore deduce 

(0|Q B (0) $ 5 fl(0)|0> = V^W), (A.7) 

as quoted in eq.(2.15). 

The renormalisation group equation for the topological susceptibility follows from the 
definition of the renormalised composite operators, eq.(2.4), and the chiral Ward identities. 
The Ward identity for the two-current Green function is 

i*? W% R (k) 4 5R (-k)\0) - 2N F (0\Q R (k) J° 5R (-k)\0) = 0. (A8) 

Combining eqs.(A.l), (A.8) and (2.4), we find straightforwardly 

(0\Q R (k) Q R (-k)\0) = Z 2 (0\Q B (k) Q B (-k)\0) + .... (A. 9) 

The dots denote the extra divergences associated with contact terms in the two-point 
Green functions of composite operators. Taking these into account (see refs.[21,7] for full 
details) we find the full RGE for X (k 2 ), 

(4 + " K)Q '^- 27 ) x(fc2) = -(^F 2/J<1,(a * )fc4 ' (A10) 

where is a new RG function. The inhomogeneous term does not contribute at zero 
momentum, however, and the required RGE (2.13) for x'(0) follows immediately. 
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Appendix B 



Decay Constants and the 77' 



We can estimate the parameter f v > appearing in the spectral expansion using the first 
Laplace QSSR, eq.(3.10). f v ' is denned by 

W% R (k)W) = iKU, (B.l) 

and is RG non-invariant. On shell (see ref.[7], Appendix D), the scale dependence is due 
entirely to the anomalous dimension 7 of the axial current so, using eq.(3.3) and expressing 
the result in terms of the QCD scale A, we may write 

Ufa) = U'^^~^ {B.2) 



where f v > is RG invariant. From the QSSR (3.10), we find the r stability starts at t c ~ 
6.5GeV 2 , while the t c stability is reached for t c larger than 9.5GeV 2 . In this region, the 
radiative corrections are about 10% of the lowest order term, while the < g 3 G 3 > one 
contributes about 10%. Under such conditions, our optimal result at r ~ 0.6GeV -2 is (see 
Fig. 6a): 

/„, ~ (24.1 ±0.6 ±3.4 ±0.3)MeV, (B.3) 

where the first error comes from < a s G 2 >, the second from A and the third from the 
range of t c values between 6.5 GeV 2 and 9.5 GeV 2 . Adding a 5% error from the unknown 
QCD terms, adding the different errors quadratically and running to the EMC scale, we 
obtain 

U\emc ^ (23.6 ± 3.5)MeV. (BA) 

This value is strongly suppressed relative to the OZI prediction of y/Gf^ for the r\' decay 
constant. 

However, as has been shown in refs.[41,7], this f v > is not the r\' decay constant measured 
in, e.g., the decay 7/ — > 77. In fact, the analogues of the current algebra formulae 

and 

UgnNN = m N gA, (B.6) 
in the flavour singlet sector are [4 1,7] 

F ^'77 + ^TT F2m2 ri' 9Grr( Q ) = -«em, (B.7) 



and 



Fg v >NN + ^F 2 ml,g GNN (Q) = 2m N G ( %\o). (B.8) 
ziv p 
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Here, F is the RG invariant decay constant denned by 



F = 



_ i 

^(/ dx i(0\T*<fi R (x) 0jj B (O)|O>) \ (5.9) 



where 0° = i^Tl^q and ((f)) = ^2(qq)- The extra terms gG-y-y an d Qgnn appearing in 
eqs.(B.7,B.8) (which are properly defined as proper vertices[41,7]) may be thought of as the 
couplings of the gluonic component of the rj'. They arise because the 7/ is not a Goldstone 
boson in the U(l) channel and so the naive current algebra extensions of eqs.(B.5,B.6) are 
not valid. At first sight, therefore, eqs.(B.7) and (B.8) are not predictive since gc-y-y and 
9gnn are unknown. However, if we follow our proposal that OZI violations are associated 
with RG non-invariant quantities we can make predictions. 

Taking eq.(B.7) first, we have shown[41] that go-y-y is RG invariant. Since in the 
OZI limit this term is absent, we therefore expect go-yy to be small, and so to a good 
approximation we predict 

4 

Fgri'-ry = ""em- (B .10) 



n 



Since F is RG invariant, we expect it to be well approximated by its OZI value V6f n . 
Experimentally (see ref.[42]), the relation (B.10) is very well satisfied. 

In eq.(B.8), on the other hand, gcNN is not RG invariant so we do not expect this 
term to be small. In fact, this equation is just a rewriting of the U(l) GT formula quoted 
in the text, for which our proposal is successful. 

An important test of our picture of the pattern of OZI breaking is therefore to evaluate 
the RG invariant decay constant F from first principles and check that it is close to the 
OZI prediction of \ / 2N F f 7r . Again, we can use QCD spectral sum rules. 

We require the zero-momentum limit $5(0) of the two-point correlation function 

d> 5 (fc 2 ) = j dxe ik -*i(0\T*4>l R {x) 0g fl (O)|O> (5.11) 

for QCD with 3 flavours and massless quarks. However, as there is a smooth behaviour 
of the two-point correlator when the common light quark mass m R goes to zero, we shall 
work (for convenience) with the RG invariant correlation function 

V 5 (k 2 ) = 4m 2 R ^ 5 (k 2 ) (B.12) 

where tur is the average of the renormalised u and d quark masses. Now, in perturbation 
theory, the difference between this flavour singlet correlation function and the correspond- 
ing non-singlet one appears only at 0(a 2 s ) from the double-triangle anomaly-type diagrams. 
Similarly for the non-perturbative condensate terms, the difference is only of 0(a 2 ) arising 
from the equivalent diagrams. Instanton-like effects appear as higher-dimension operators. 
So, at the order we are working, we can simply use the expression for the isotriplet (pion) 
correlation function in QCD discussed in the literature [20]. 
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The first Laplace sum rule to two-loops reads [20]: 



Jo 



1 



dte- tT -lm^ 5 (t) ~ , 

TV ZTT 



^fm 2 (r)|r- 2 (l-e—(l + t c r)) 



2/11 2 / „ #A 7!/? 2 

1 - ALU +27£ "Al 72 - 7i Aj +2 ^r logL 



+ 



7T 2 896 o , . o 

-(a s G ) + — 7T pa s {uu) r 



(5.13) 



where L = — log r A 2 and [20] : 

po^-u-u) 2 ~ (3.8±2.0)10" 4 GeV 6 

m(r) = J(m u + m d )(r) ~ (-^ log rA 2 )^ l//Sl (12.1 ± 1.0)MeV (S.14) 

— ■ — ■ 



As before, we parametrise the spectral function keeping only the lowest (77') resonance, 



i.e. 



1 



TV 



Im# 5 (£) = 2mtf 2 5(t - fh 2 ,) + "QCD continuum" 9(t - t c ), 



where the unknown parameter /, which is defined by 

277^(01^17/) = V2fm 2 v , : 



(B.15) 



(5.16) 



can be estimated from the sum rule (B.13). We study the r and t c behaviours of / in Fig. 7. 
The r stability starts for t c ~ 4GeV 2 , while stability in t c appears above t c ~ 7GeV 2 , a 
range which is equal to the one for the correlation function for Q(x). The value for the r 
stability of about 0.9GeV -2 is typical of light quark correlation functions. At the minimum, 
we obtain 

/ = v /iV> (5.55 ±0.08 ±0.65 ±0.35 ±0.06 ±0.03) MeV, (5.17) 

where the errors come respectively from t c , A, 771, (a s G 2 ) and pa s {uu) 2 . Adding these 
errors quadratically, we deduce 



/ = 0VV (5.55 ± 0.75)MeV. 



(5.18) 



With this value for /, we are now able to estimate ^5(0) itself using a second Laplace 
sum rule [43, 25]: 



* 5 (0) 



dt 
— < 
t 



m 2 (r) 



jr-^l-e-*^) 



+ 



11 9 2/ „ fo\ 7x/? 2 

p 1 LKY + 2 ^-W 1 V 2 - 11 Wj +2 W losL 



tv 2 1 896 o 2 
-(a s G ) + - — tv pa s (uu) r 



(5.19) 
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where 71, 72 are the coefficients in the anomalous dimension for the light quark mass. For 
three flavours, 71 = 2 and 72 = 91/12. The sum rule analysis of this quantity shows a 
strong t c dependence and the r stability only appears at unrealistic values of t c larger than 
8GeV 2 . In order to circumvent this difficulty, we work with a combination of the sum rules 
(B.13) and (B.19) which has been used successfully in the past for measuring the deviation 
from pion and kaon PC AC to a good accuracy [43]. The combined sum rule reads: 



* 5 (0) ^ 



dt 



3N F 



TO' 



+ T 



n 2/ ~ 9 7i/? 2l r 

ALVT +27£ -Al 72 ~ 7i Aj +2 ^r logL 



27T. ~. 3 896 3 ._ ' 
y (a s G ) + - —n pa s (uu) t 



(5.20) 



This sum rule is studied in Fig. 8. The position of the stability is almost insensitive to 
the value of t c due to some cancellations amongst the perturbative terms. However, this 
feature also implies that the stability is obtained at values of r larger than in the previous 
cases making the result sensitive to the errors on the four-quark condensates, which affects 
the accuracy of the result. We deduce 



* 5 (0) ~ N F (3.70 ± 0.90 ± 0.30 ± 0.70 ± 2.00) 10" 6 GeV 4 , 



(5.21) 



where the errors are due to /, A, (a s G 2 ) and pa s (uu) 2 . Adding these errors quadratically, 

we obtain 

y/V 5 (0) ^ ^/Nf (1.92±0.53)10- 3 GeV 2 . (5.22) 

Using this value in eq.(B.9) (with rh v >), after multiplying the numerator and denominator 
by the overall 2m r factor and using Dashen's formula for iur^r), we finally find 



F ~ (1.55 ± 0.43) v / 2iV> /*■, 



(5.23) 



to be compared with the OZI prediction of \f2Np f n . 

This result is again in broad agreement with our expectations, although of course the 
errors are much too large to draw a definitive conclusion. Nevertheless, this confirmation 
can be taken as providing extra support for the reliability of the estimate in the text for 
X'(0). 
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Figure captions 



Fig. 1 The two-current matrix element (-ZV|J M (g) J u (—q)\N). 

Fig. 2 Decomposition of the matrix element into a composite operator propagator (de- 
noted by the double line) and a proper vertex (hatched). 

Fig. 3 The original parton model representation of the scattering amplitude. 
Fig. 4 The QCD-improved parton model representation. 

Fig. 5 a) r-behaviour of a/ x'(0) for different values of the continuum threshold t c . 
b) behaviour of different r-minima versus t c . 

Fig. 6 a) As Fig. 5a for the parameter f Vl . 

b) FESR prediction of vV(0) versus t c for different values of f v >. 

Fig. 7 As Fig. 5a for the parameter /. 

Fig. 8 As Fig. 5a for # 5 (0). 
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